FixNCut: single-cell genomics through reversible tissue fixation and dissociation.

The use of single-cell technologies for clinical applications requires disconnecting sampling from downstream processing steps. Early sample preservation can further increase robustness and reproducibility by avoiding artifacts introduced during specimen handling. We present FixNCut, a methodology for the reversible fixation of tissue followed by dissociation that overcomes current limitations. We applied FixNCut to human and mouse tissues to demonstrate the preservation of RNA integrity, sequencing library complexity, and cellular composition, while diminishing stress-related artifacts. Besides single-cell RNA sequencing, FixNCut is compatible with multiple single-cell and spatial technologies, making it a versatile tool for robust and flexible study designs.

Multi-organ single-cell transcriptomics of immune cells uncovered organ-specific gene expression and functions.

Despite the wealth of publicly available single-cell datasets, our understanding of distinct resident immune cells and their unique features in diverse human organs remains limited. To address this, we compiled a meta-analysis dataset of 114,275 CD45+ immune cells sourced from 14 organs in healthy donors. While the transcriptome of immune cells remains relatively consistent across organs, our analysis has unveiled organ-specific gene expression differences (GTPX3 in kidney, DNTT and ACVR2B in thymus). These alterations are linked to different transcriptional factor activities and pathways including metabolism. TNF-α signaling through the NFkB pathway was found in several organs and immune compartments. The presence of distinct expression profiles for NFkB family genes and their target genes, including cytokines, underscores their pivotal role in cell positioning. Taken together, immune cells serve a dual role: safeguarding the organs and dynamically adjusting to the intricacies of the host organ environment, thereby actively contributing to its functionality and overall homeostasis.

An atlas of cells in the human tonsil.

Palatine tonsils are secondary lymphoid organs (SLOs) representing the first line of immunological defense against inhaled or ingested pathogens. We generated an atlas of the human tonsil composed of >556,000 cells profiled across five different data modalities, including single-cell transcriptome, epigenome, proteome, and immune repertoire sequencing, as well as spatial transcriptomics. This census identified 121 cell types and states, defined developmental trajectories, and enabled an understanding of the functional units of the tonsil. Exemplarily, we stratified myeloid slan-like subtypes, established a BCL6 enhancer as locally active in follicle-associated T and B cells, and identified SIX5 as putative transcriptional regulator of plasma cell maturation. Analyses of a validation cohort confirmed the presence, annotation, and markers of tonsillar cell types and provided evidence of age-related compositional shifts. We demonstrate the value of this resource by annotating cells from B cell-derived mantle cell lymphomas, linking transcriptional heterogeneity to normal B cell differentiation states of the human tonsil.

Convert-Pheno: A software toolkit for the interconversion of standard data models for phenotypic data.

Efficient sharing and integration of phenotypic data is crucial for advancing biomedical research and enhancing patient outcomes in precision medicine and public health. To achieve this, the health data community has developed standards to promote the harmonization of variable names and values. However, the use of diverse standards across different research centers can hinder progress. Here we present Convert-Pheno, an open-source software toolkit that enables the interconversion of common data models for phenotypic data such as Beacon v2 Models, CDISC-ODM, OMOP-CDM, Phenopackets v2, and REDCap. Along with the software, we have created a detailed documentation that includes information on deployment and installation.

Development of a novel streamlined workflow (AACRE) and database (inCREDBle) for genomic analysis of carbapenem-resistant Enterobacterales.

In response to the threat of increasing antimicrobial resistance, we must increase the amount of available high-quality genomic data gathered on antibiotic-resistant bacteria. To this end, we developed an integrated pipeline for high-throughput long-read sequencing, assembly, annotation and analysis of bacterial isolates and used it to generate a large genomic data set of carbapenemase-producing Enterobacterales (CPE) isolates collected in Spain. The set of 461 isolates were sequenced with a combination of both Illumina and Oxford Nanopore Technologies (ONT) DNA sequencing technologies in order to provide genomic context for chromosomal loci and, most importantly, structural resolution of plasmids, important determinants for transmission of antimicrobial resistance. We developed an informatics pipeline called Assembly and Annotation of Carbapenem-Resistant Enterobacteriaceae (AACRE) for the full assembly and annotation of the bacterial genomes and their complement of plasmids. To explore the resulting genomic data set, we developed a new database called inCREDBle that not only stores the genomic data, but provides unique ways to filter and compare data, enabling comparative genomic analyses at the level of chromosomes, plasmids and individual genes. We identified a new sequence type, ST5000, and discovered a genomic locus unique to ST15 that may be linked to its increased spread in the population. In addition to our major objective of generating a large regional data set, we took the opportunity to compare the effects of sample quality and sequencing methods, including R9 versus R10 nanopore chemistry, on genome assembly and annotation quality. We conclude that converting short-read and hybrid microbial sequencing and assembly workflows to the latest nanopore chemistry will further reduce processing time and cost, truly enabling the routine monitoring of resistance transmission patterns at the resolution of complete chromosomes and plasmids.

Genetic newborn screening and digital technologies: A project protocol based on a dual approach to shorten the rare diseases diagnostic path in Europe.

Since 72% of rare diseases are genetic in origin and mostly paediatrics, genetic newborn screening represents a diagnostic "window of opportunity". Therefore, many gNBS initiatives started in different European countries. Screen4Care is a research project, which resulted of a joint effort between the European Union Commission and the European Federation of Pharmaceutical Industries and Associations. It focuses on genetic newborn screening and artificial intelligence-based tools which will be applied to a large European population of about 25.000 infants. The neonatal screening strategy will be based on targeted sequencing, while whole genome sequencing will be offered to all enrolled infants who may show early symptoms but have resulted negative at the targeted sequencing-based newborn screening. We will leverage artificial intelligence-based algorithms to identify patients using Electronic Health Records (EHR) and to build a repository "symptom checkers" for patients and healthcare providers. S4C will design an equitable, ethical, and sustainable framework for genetic newborn screening and new digital tools, corroborated by a large workout where legal, ethical, and social complexities will be addressed with the intent of making the framework highly and flexibly translatable into the diverse European health systems.

A chromosome-level reference genome for the common octopus, Octopus vulgaris (Cuvier, 1797).

Cephalopods are emerging animal models and include iconic species for studying the link between genomic innovations and physiological and behavioral complexities. Coleoid cephalopods possess the largest nervous system among invertebrates, both for cell counts and brain-to-body ratio. Octopus vulgaris has been at the center of a long-standing tradition of research into diverse aspects of cephalopod biology, including behavioral and neural plasticity, learning and memory recall, regeneration, and sophisticated cognition. However, no chromosome-scale genome assembly was available for O. vulgaris to aid in functional studies. To fill this gap, we sequenced and assembled a chromosome-scale genome of the common octopus, O. vulgaris. The final assembly spans 2.8 billion basepairs, 99.34% of which are in 30 chromosome-scale scaffolds. Hi-C heatmaps support a karyotype of 1n = 30 chromosomes. Comparisons with other octopus species' genomes show a conserved octopus karyotype and a pattern of local genome rearrangements between species. This new chromosome-scale genome of O. vulgaris will further facilitate research in all aspects of cephalopod biology, including various forms of plasticity and the neural machinery underlying sophisticated cognition, as well as an understanding of cephalopod evolution.

Germ line variants in patients with acute myeloid leukemia without a suspicion of hereditary hematologic malignancy syndrome.

Germ line predisposition in acute myeloid leukemia (AML) has gained attention in recent years because of a nonnegligible frequency and an impact on management of patients and their relatives. Risk alleles for AML development may be present in patients without a clinical suspicion of hereditary hematologic malignancy syndrome. In this study we investigated the presence of germ line variants (GVs) in 288 genes related to cancer predisposition in 47 patients with available paired, tumor-normal material, namely bone marrow stroma cells (n = 29), postremission bone marrow (n = 17), and saliva (n = 1). These patients correspond to 2 broad AML categories with heterogeneous genetic background (AML myelodysplasia related and AML defined by differentiation) and none of them had phenotypic abnormalities, previous history of cytopenia, or strong cancer aggregation. We found 11 pathogenic or likely pathogenic variants, 6 affecting genes related to autosomal dominant cancer predisposition syndromes (ATM, DDX41, and CHEK2) and 5 related to autosomal recessive bone marrow failure syndromes (FANCA, FANCM, SBDS, DNAJC21, and CSF3R). We did not find differences in clinical characteristics nor outcome between carriers of GVs vs noncarriers. Further studies in unselected AML cohorts are needed to determine GV incidence and penetrance and, in particular, to clarify the role of ATM nonsense mutations in AML predisposition.

Dissecting the Origin of Heterogeneity in Uterine and Ovarian Carcinosarcomas.

Gynecologic carcinosarcomas (CS) are biphasic neoplasms composed of carcinomatous (C) and sarcomatous (S) malignant components. Because of their rarity and histologic complexity, genetic and functional studies on CS are scarce and the mechanisms of initiation and development remain largely unknown. Whole-genome analysis of the C and S components reveals shared genomic alterations, thus emphasizing the clonal evolution of CS. Reconstructions of the evolutionary history of each tumor further reveal that C and S samples are composed of both ancestral cell populations and component-specific subclones, supporting a common origin followed by distinct evolutionary trajectories. However, while we do not find any recurrent genomic features associated with phenotypic divergence, transcriptomic and methylome analyses identify a common mechanism across the cohort, the epithelial-to-mesenchymal transition (EMT), suggesting a role for nongenetic factors in inflicting changes to cellular fate. Altogether, these data accredit the hypothesis that CS tumors are driven by both clonal evolution and transcriptomic reprogramming, essential for susceptibility to transdifferentiation upon encountering environmental cues, thus linking CS heterogeneity to genetic, transcriptomic, and epigenetic influences. Significance: We have provided a detailed characterization of the genomic landscape of CS and identified EMT as a common mechanism associated with phenotypic divergence, linking CS heterogeneity to genetic, transcriptomic, and epigenetic influences.

A critical spotlight on the paradigms of FFPE-DNA sequencing.

In the late 19th century, formalin fixation with paraffin-embedding (FFPE) of tissues was developed as a fixation and conservation method and is still used to this day in routine clinical and pathological practice. The implementation of state-of-the-art nucleic acid sequencing technologies has sparked much interest for using historical FFPE samples stored in biobanks as they hold promise in extracting new information from these valuable samples. However, formalin fixation chemically modifies DNA, which potentially leads to incorrect sequences or misinterpretations in downstream processing and data analysis. Many publications have concentrated on one type of DNA damage, but few have addressed the complete spectrum of FFPE-DNA damage. Here, we review mitigation strategies in (I) pre-analytical sample quality control, (II) DNA repair treatments, (III) analytical sample preparation and (IV) bioinformatic analysis of FFPE-DNA. We then provide recommendations that are tested and illustrated with DNA from 13-year-old liver specimens, one FFPE preserved and one fresh frozen, applying target-enriched sequencing. Thus, we show how DNA damage can be compensated, even when using low quantities (50 ng) of fragmented FFPE-DNA (DNA integrity number 2.0) that cannot be amplified well (Q129 bp/Q41 bp = 5%). Finally, we provide a checklist called 'ERROR-FFPE-DNA' that summarises recommendations for the minimal information in publications required for assessing fitness-for-purpose and inter-study comparison when using FFPE samples.

Remote visualization of large-scale genomic alignments for collaborative clinical research and diagnosis of rare diseases.

The Solve-RD project objectives include solving undiagnosed rare diseases (RD) through collaborative research on shared genome-phenome datasets. The RD-Connect Genome-Phenome Analysis Platform (GPAP), for data collation and analysis, and the European Genome-Phenome Archive (EGA), for file storage, are two key components of the Solve-RD infrastructure. Clinical researchers can identify candidate genetic variants within the RD-Connect GPAP and, thanks to the developments presented here as part of joint ELIXIR activities, are able to remotely visualize the corresponding alignments stored at the EGA. The Global Alliance for Genomics and Health (GA4GH) htsget streaming application programming interface (API) is used to retrieve alignment slices, which are rendered by an integrated genome viewer (IGV) instance embedded in the GPAP. As a result, it is no longer necessary for over 11,000 datasets to download large alignment files to visualize them locally. This work highlights the advantages, from both the user and infrastructure perspectives, of implementing interoperability standards for establishing federated genomics data networks.

Machine learning and multi-omics data in chronic lymphocytic leukemia: the future of precision medicine?

Chronic lymphocytic leukemia is a complex and heterogeneous hematological malignancy. The advance of high-throughput multi-omics technologies has significantly influenced chronic lymphocytic leukemia research and paved the way for precision medicine approaches. In this review, we explore the role of machine learning in the analysis of multi-omics data in this hematological malignancy. We discuss recent literature on different machine learning models applied to single omic studies in chronic lymphocytic leukemia, with a special focus on the potential contributions to precision medicine. Finally, we highlight the recently published machine learning applications in multi-omics data in this area of research as well as their potential and limitations.

DNA methylation in human gastric epithelial cells defines regional identity without restricting lineage plasticity.

BACKGROUND: Epigenetic modifications in mammalian DNA are commonly manifested by DNA methylation. In the stomach, altered DNA methylation patterns have been observed following chronic Helicobacter pylori infections and in gastric cancer. In the context of epigenetic regulation, the regional nature of the stomach has been rarely considered in detail. RESULTS: Here, we establish gastric mucosa derived primary cell cultures as a reliable source of native human epithelium. We describe the DNA methylation landscape across the phenotypically different regions of the healthy human stomach, i.e., antrum, corpus, fundus together with the corresponding transcriptomes. We show that stable regional DNA methylation differences translate to a limited extent into regulation of the transcriptomic phenotype, indicating a largely permissive epigenetic regulation. We identify a small number of transcription factors with novel region-specific activity and likely epigenetic impact in the stomach, including GATA4, IRX5, IRX2, PDX1 and CDX2. Detailed analysis of the Wnt pathway reveals differential regulation along the craniocaudal axis, which involves non-canonical Wnt signaling in determining cell fate in the proximal stomach. By extending our analysis to pre-neoplastic lesions and gastric cancers, we conclude that epigenetic dysregulation characterizes intestinal metaplasia as a founding basis for functional changes in gastric cancer. We present insights into the dynamics of DNA methylation across anatomical regions of the healthy stomach and patterns of its change in disease. Finally, our study provides a well-defined resource of regional stomach transcription and epigenetics.

Tackling the translational challenges of multi-omics research in the realm of European personalised medicine: A workshop report.

Personalised medicine (PM) presents a great opportunity to improve the future of individualised healthcare. Recent advances in -omics technologies have led to unprecedented efforts characterising the biology and molecular mechanisms that underlie the development and progression of a wide array of complex human diseases, supporting further development of PM. This article reflects the outcome of the 2021 EATRIS-Plus Multi-omics Stakeholder Group workshop organised to 1) outline a global overview of common promises and challenges that key European stakeholders are facing in the field of multi-omics research, 2) assess the potential of new technologies, such as artificial intelligence (AI), and 3) establish an initial dialogue between key initiatives in this space. Our focus is on the alignment of agendas of European initiatives in multi-omics research and the centrality of patients in designing solutions that have the potential to advance PM in long-term healthcare strategies.

A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome.

Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research.

The ulcerative colitis-associated gene FUT8 regulates the quantity and quality of secreted mucins.

Mucins are the main macrocomponents of the mucus layer that protects the digestive tract from pathogens. Fucosylation of mucins increases mucus viscoelasticity and its resistance to shear stress. These properties are altered in patients with ulcerative colitis (UC), which is marked by a chronic inflammation of the distal part of the colon. Here, we show that levels of Fucosyltransferase 8 (FUT8) and specific mucins are increased in the distal inflamed colon of UC patients. Recapitulating this FUT8 overexpression in mucin-producing HT29-18N2 colonic cell line increases delivery of MUC1 to the plasma membrane and extracellular release of MUC2 and MUC5AC. Mucins secreted by FUT8 overexpressing cells are more resistant to removal from the cell surface than mucins secreted by FUT8-depleted cells (FUT8 KD). FUT8 KD causes intracellular accumulation of MUC1 and alters the ratio of secreted MUC2 to MUC5AC. These data fit well with the Fut8-/- mice phenotype, which are protected from UC. Fut8-/- mice exhibit a thinner proximal colon mucus layer with an altered ratio of neutral to acidic mucins. Together, our data reveal that FUT8 modifies the biophysical properties of mucus by controlling levels of cell surface MUC1 and quantity and quality of secreted MUC2 and MUC5AC. We suggest that these changes in mucus viscoelasticity likely facilitate bacterial-epithelial interactions leading to inflammation and UC progression.

Detection of early seeding of Richter transformation in chronic lymphocytic leukemia.

Richter transformation (RT) is a paradigmatic evolution of chronic lymphocytic leukemia (CLL) into a very aggressive large B cell lymphoma conferring a dismal prognosis. The mechanisms driving RT remain largely unknown. We characterized the whole genome, epigenome and transcriptome, combined with single-cell DNA/RNA-sequencing analyses and functional experiments, of 19 cases of CLL developing RT. Studying 54 longitudinal samples covering up to 19 years of disease course, we uncovered minute subclones carrying genomic, immunogenetic and transcriptomic features of RT cells already at CLL diagnosis, which were dormant for up to 19 years before transformation. We also identified new driver alterations, discovered a new mutational signature (SBS-RT), recognized an oxidative phosphorylation (OXPHOS)high-B cell receptor (BCR)low-signaling transcriptional axis in RT and showed that OXPHOS inhibition reduces the proliferation of RT cells. These findings demonstrate the early seeding of subclones driving advanced stages of cancer evolution and uncover potential therapeutic targets for RT.

The RD-Connect Genome-Phenome Analysis Platform: Accelerating diagnosis, research, and gene discovery for rare diseases.

Rare disease patients are more likely to receive a rapid molecular diagnosis nowadays thanks to the wide adoption of next-generation sequencing. However, many cases remain undiagnosed even after exome or genome analysis, because the methods used missed the molecular cause in a known gene, or a novel causative gene could not be identified and/or confirmed. To address these challenges, the RD-Connect Genome-Phenome Analysis Platform (GPAP) facilitates the collation, discovery, sharing, and analysis of standardized genome-phenome data within a collaborative environment. Authorized clinicians and researchers submit pseudonymised phenotypic profiles encoded using the Human Phenotype Ontology, and raw genomic data which is processed through a standardized pipeline. After an optional embargo period, the data are shared with other platform users, with the objective that similar cases in the system and queries from peers may help diagnose the case. Additionally, the platform enables bidirectional discovery of similar cases in other databases from the Matchmaker Exchange network. To facilitate genome-phenome analysis and interpretation by clinical researchers, the RD-Connect GPAP provides a powerful user-friendly interface and leverages tens of information sources. As a result, the resource has already helped diagnose hundreds of rare disease patients and discover new disease causing genes.

Emergence of 16S rRNA methyltransferases among carbapenemase-producing Enterobacterales in Spain studied by whole-genome sequencing.

The emergence of 16S rRNA methyltransferases (RMTs) in Gram-negative pathogens bearing other clinically relevant resistance mechanisms, such as carbapenemase-producing Enterobacterales (CPE), is becoming an alarming concern. We investigated the prevalence, antimicrobial susceptibility, resistance mechanisms, molecular epidemiology and genetic support of RMTs in CPE isolates from Spain. This study included a collection of 468 CPE isolates recovered during 2018 from 32 participating Spanish hospitals. MICs were determined using the broth microdilution method, the agar dilution method (fosfomycin) or MIC gradient strips (plazomicin). All isolates were subjected to hybrid whole-genome sequencing (WGS). Sequence types (STs), core genome phylogenetic relatedness, horizontally acquired resistance mechanisms, plasmid analysis and the genetic environment of RMTs were determined in silico from WGS data in all RMT-positive isolates. Among the 468 CPE isolates evaluated, 24 isolates (5.1%) recovered from nine different hospitals spanning five Spanish regions showed resistance to all aminoglycosides and were positive for an RMT (21 RmtF, 2 ArmA and 1 RmtC). All RMT-producers showed high-level resistance to all aminoglycosides, including plazomicin, and in most cases exhibited an extensively drug-resistant susceptibility profile. The RMT-positive isolates showed low genetic diversity and were global clones of Klebsiella pneumoniae (ST147, ST101, ST395) and Enterobacter cloacae (ST93) bearing blaOXA-48, blaNDM-1 or blaVIM-1 carbapenemase genes. RMTs were harboured in five different multidrug resistance plasmids and linked to efficient mobile genetic elements. Our findings highlight that RMTs are emerging among clinical CPE isolates from Spain and their spread should be monitored to preserve the future clinical utility of aminoglycosides and plazomicin.